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Statistical ensembles of flexible two-dimensional fluid membranes arise naturally in the 
description of many physical systems. Typically one encounters such systems in a regime 
of low tension but high stiffness against bending, which is just the opposite of the regime 
described by the Polyakov string. We study a class of couplings between membrane shape 
and in-plane order which break 3-space parity invariance. Remarkably there is only one 
such allowed coupling (up to boundary terms) ; this term will be present for any lipid bilayer 
composed of tilted chiral molecules. We calculate the renormalizat ion-group behavior of 
this relevant coupling in a simplified model and show how thermal fluctuations effectively 
reduce it in the infrared. 
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Statistical ensembles of random geometrical shapes pervade theoretical physics. Ini- 
tially one-dimensional curves in space were most thoroughly studied due to their ease of 
description and the many applications of such ensembles to the conformation and dynamics 
of polymers, but today ensembles of two-dimensional membranes are at least as impor- 
tant. Physical realizations of such surfaces include lipid bilayers and surfactant films, which 
spontaneously self-assemble from amphiphilic molecules in solution or at fluid interfaces 
(for reviews see [J]). More speculative applications as diverse as the 3d Ising model and 
other 3d phase transitions, cosmic strings, flux tubes in QCD, and models of elementary 
particles all rest upon the key property that the important physical degrees of freedom in 
these problems are shapes with no preferred choice of coordinates. The condition that 
coordinate choice be immaterial greatly constrains the possible forms of the statistical 
weights in these systems, leading to very few independent couplings and hence physically 
simple models. 

In this letter we will study a model appropriate for the description of tilted lipid 
bilayers (e.g. the lamellar L^* phases of lyotropics or S c * phases of smectics), though we 
think the analysis is potentially interesting in other contexts as well. Bilayers are typically 
rigid, that is their resistance to bending is characterized by a dimensionless quantity R = 
K /k B T (see below) which is greater than one. Accordingly, we will carry out a perturbative 
expansion about R — > oo, the high-stiffness, low-temperature limit. Bilayers with free 
boundary conditions also typically adjust themselves to zero effective surface tension 0, 
and we will also work in this limit. II 

We can summarize our logic as follows (further details will appear elsewhere). At 
length scales far longer than the size of the constituent molecules a continuum description 
becomes appropriate. In our nearly flat regime the important degrees of freedom are the 
elastic (or "Goldstone" ) modes corresponding to transverse undulations as well as director 

1 Sheets of molecules frozen into fixed lattices form "tethered" membranes which do have a 
preferred coordinate system, much like elastic solids. Here we study only the opposite case ( "fluid" 
membranes) . 

2 This is opposite to the usual Polyakov string, in which tension dominates and Rq is effectively 
zero. Refs. ||][|| have shown in a related model (hexatic membrane) that with in-plane order 
the stiffness can stabilize at a large value without running to zero at long scales. In any case 
the high-stiffness limit is appropriate for a system viewed on scales shorter than its persistence 
length, if any. In this paper, we will not study the running of the various stiffness couplings at 
all, concentrating instead on the behavior of chirality. 
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fluctuations if in-plane order develops. We will find only a few allowed couplings of these 
modes, as we expect in any elastic system at very long wavelengths: the system forgets 
most of the details about its constituents. Somewhat surprisingly, however, the membranes 
with in-plane order can remember at long scales whether or not their constituent molecules 
are chiral, even though this chirality is often a rather subtle property of the amphiphiles. 
Gross chiral behavior has long been seen in monolayers, where one gets pinwheel domains 
H ,! as well as in flexible membranes which can form helical ribbons || whose sense depends 
on the constituent molecules' handedness 0. We will see that this memory follows from 
the existence of an allowed parity- violating term in the free energy, which couples in-plane 
order to shape. |§ When in-plane order is thermally destroyed the theory admits no such 
terms; the system cannot express the chirality of its constituents even if present. Indeed, 
experimentally a loss of gross chiral structure does seem to accompany the chain-melting 
transition, and chiral structures do not form at all above this temperature || J7| . 

Thermal fluctuations are often important in membranes at room temperature where 
R < 40 is not too close to infinity. It is well known that without in-plane order such 
fluctuations induce a logarithmically scale-dependent softening of the effective stiffness R 



in the infrared |T(J. What we will show is that the unique bulk chiral coupling, if present, 
similarly suffers a logarithmic renormalization The RG behavior of this chiral term seems 
not to have been studied before. While some of our analysis will reproduce others' results, 
we hope that our unified treatment of allowed free energy terms will clarify some of the 
important symmetries; we have also tried to clear up several subtleties in the fluctuation 



problem |TT| . Finally, the renormalization of chirality will affect the average shapes taken 
on by membranes, since the chiral coupling helps determine those shapes. We will derive 
an anomalous scaling relation for the radius of helical ribbons as chirality varies which 
departs from the mean-field formula |E| and may be experimentally testable.!! 

We begin by considering in greater detail the elastic modes of our system. A membrane 
made of molecules which slip around each other in some average 2d locus with no in- 
plane order at all can be described just by specifying a mathematical 2-surface in 
3-space. We can choose coordinates £\ i = 1,2 for our surface, but we must remember 



3 We must be careful not to confuse the explicit breaking of parity symmetry due to chiral 
molecules with the spontaneous breaking of parity symmetry (R. Bruinsma and J. Selinger, to 
appear) seen in pinwheel domains made from achiral molecules (X. Qiu et al, Phys. Rev. Lett. 
67 (1991) 703). 

We thank J. Toner for suggesting this consequence. 
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that this choice is arbitrary. The choice of a flat reference surface breaks one transverse 
translation symmetry, giving one "undulation" mode. We can visualize this mode as height 
fluctuations from a horizontal plane. At high temperatures this is the only soft mode 
we expect. At lower temperatures in-plane order can develop, leading to additional soft 
modes. While experiments do not seem to see full crystalline in-plane order in (hydrated, 
unpolymerized) membranes [|| , and such order is disfavored on theoretical grounds , still 
orientational quasi-long-range order can survive, giving rise to an angular elastic mode. 
We will consider only temperatures well within the ordered phase. Hence we can visualize 
the corresponding order parameter as a unit vector field m(£) tangent to our 2-surface. 

For hexatic in-plane order m is defined only up to rotations by 2n/6. When the 
constituent molecules are tilted from the surface normal they break rotation invariance 
completely, so m has no periodic identifications, but a new subtlety arises instead. To 
define given a configuration of molecules, we must project the molecule axis down 

to the midplane of the membrane. Since the two sides of the bilayers are equivalent, the 
overall sign of is not fixed until we choose one side, or equivalently, one of the two 
normal vectors n(£); with this choice we may use the convention where the projection of 
the molecule on the outward-facing layer defines m. Since there is no preferred choice of 
normal, each term of our free energy must be unchanged under the substitution rh i— > — m, 
ft I— > — n.i 

Our surface gets a few standard tensor fields. It inherits a metric g- = d { x ■ djX, 
a corresponding covariant derivative V, and a volume form d 2 £y/g. Once we choose a 
normal n, we also get a second fundamental form = n ■ V \djX and an alternating 

tensor e y = e lja n a = y/g I ^ M . Here e ija = d^d^e^ where e bca is the usual 

\ / ij 

alternating symbol, e 123 = +1. Note that and change sign under the change of 
normal n h- > — h, m h-> — m discussed above. 

Let us enumerate all allowed free energy terms. We assume all nonlocal interactions 
are absent or screened to a length scale shorter than the scale of interest. Then we simply 
seek all local terms in x{£) and m(£) with the above symmetries, relevant or marginal in 
the low-temperature or high-stiffness expansion. We will not seek to derive such terms 
from some 3-dimensional liquid crystal free energy, but instead simply construct them 
directly from the ingredients listed above. The relevance of an operator about the weakly 

5 For monolayers or closed vesicles there is a preferred choice of normal and we should not 
impose this symmetry. 
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fluctuating R — > oo fixed point depends on its naive engineering dimension. The only 
short-distance cutoff in the problem is the characteristic size A -1 of the constituents. A -1 
is a 3-space length, not a parameter-space length as in ordinary 2d quantum field theory, 
so when counting the dimensions of operators we should consider only their behavior under 
rescaling x (not h and m are unit vectors and hence dimensionless. 

We can now list all the independent relevant and marginal terms allowed in the free 
energy: 



#1 



Kr 



& z i^fg {K\Y + 7l (V J m J )(V l m ,) + 7 2 (m • Vm l )(m ■ VmJ 



(1) 



H 2 =^ J d 2 ^ [a-yfh ■K-K-m + a 2 (m ■ K ■ m)(K\) 

+ a 3 (m ■ K ■ m) 2 + p^rh ■ K ■ m)(V • m) + f3 2 (K\)V ■ m 
+f3 3 K l j \7 l m j + f3 4 m e K ei m J V jm l ] , 



(2) 
(3) 



In these formulas we raise and lower i,j indices using g i -. Most of these terms have already 
been discussed by Helfrich and Prost JT^]; see also [13]. H 1 is the usual Canham-Helfrich 
elastic energy [|n[ with zero tension||, plus the covariant form of an X-Y model energy. 
H 2 contains various anisotropies in the bending energy due to the tilt of the constituents. 
r y i ,a i ,(3 i are dimensionless numbers we take to be 0(1), while k is an energy scale we 
take much larger than the temperature T (we set Boltzmann's constant = 1). Every term 
of H x 2 (every nonchiral term) is marginal. is the only allowed bulk chiral term; it 
is relevant. As mentioned earlier, there are no available relevant or marginal bulk chiral 
terms involving only shape (no tilt). 

In the enumeration (0)^(0) we have dropped all total derivatives, including for instance 
the Gaussian curvature and the covariant form of the V x rh term [^5|. Such terms are 
important for systems with boundaries (e.g. ribbons) or defects (e.g. rippled phases), but 
they will not affect our calculation of the renormalization of Cq. (Near defects one should 
also allow rh to vary in magnitude ||16||.) What is remarkable about (III) — (H) is that while 



6 One can reproduce this conclusion by carrying out the traditional analysis of divergences of 
Feynman diagrams. This power counting and the nature of the corresponding continuum limit 
differ from the procedure used in the Polyakov string where e.g. x a (£) are considered as three 
scalar fields and hence dimensionless. 
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there are a number of allowed couplings, still the list is rather short, especially compared 
to a 2d scalar at its trivial fixed point, for which infinitely many marginal terms exist. 

We will consider a limit in which R = k /T is large but Cq/T is very small compared 
to the cutoff. This is reasonable since typically chirality is a very minor feature of the long 
chain amphiphiles. Thus, we will carry all our calculations out to lowest nontrivial order 
in the loop-counting parameter (k ) _1 and the chiral coupling Cq. 

We will calculate the effects of thermal fluctuations in saddle-point (one-loop) approx- 
imation. Before we can do the required integrals, however, we must fix certain coordinate- 
choice redundancies in our description of a surface. Our problem is that we have five 
variables £?(£), rh(^) in fll])-©, while as we have seen only two are truly independent. We 
must first write x in terms of one independent field it(£). For our purposes, the easiest 
choice is "Monge gauge:" = £ 2 , «(£)), since nearly-flat surfaces can be expanded 
easily in powers of u. Next we may write m as m(£) = e*i(£) cos#(£) +e 2 (£) sin#(£), where 
{e l5 e 2 } are a field of orthonormal tangent vectors to the surface. e Q depend on #(£), but x 
does not fully determine them: we must fix an 0(2) gauge freedom. Monge gauge has the 
pleasant feature that we may choose e a = ejc^x, where e Q l = 8 l a — \d(u,d a u-\- 0(u 3 ). This 
expression is not covariant because Monge gauge is not. From now on we will raise and 
lower indices using the flat metric 8^ ; all factors will be shown explicitly. Similarly, we 
convert index type using 8 l a ; all e a z factors will be shown explicitly. The fact that we can 
choose a frame with no 0{u) terms will make Monge gauge very convenient. 

We now have all the necessary ingredients. When viewing the system on a scale 
L ^> A -1 we may forget about irrelevant couplings; moreover all the effects of fluctuations 
on scales between A -1 and (6A) _1 > A -1 may be summarized by readjusting the values 
of our couplings, since we took care to include all allowed terms and our cutoff respects 
the symmetries. Our strategy is to expand u = u + h, 9 = 9 + where u, 9 have only 
long- wavelength components, while h, ( have only wavenumbers greater than bA. @ We 
expand the free energy H in h, C about u, 9, find the quadratic terms in h, C, and integrate 

7 Since (6A) _1 is still a very short scale we are allowed to use spin- wave approximation for 
9, even though on long scales an X-Y model does not have true long-range order. We leave to 
the renormalization group the task of summing the infinite string of logarithms giving the true 
long-range behavior of 6 (cf. |fT^| ). 
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the fast modes of h, £ in Gaussian approximation to get the one-loop effective "action" 
H eS [u, 9] (see e.g. |T^]).§ Finally we read off the renormalized chiral coupling c* ff . 

To keep our formulas manageable we will truncate our model (0)-©, retaining only 
the first two terms of H 1 and of course (for more details see We expect that the 

isotropic terms retained will give a good qualitative guide to the effects of fluctuations on 
chirality. The terms we have omitted will not be generated to the order we are working, 
since they lack an extra 0(2) symmetry of the first two terms; nor can the chiral term 
induce them to first order in Cq. 

We want to pick off from H eS the renormalized coefficient of \ J y/gfh^e^ ' K^m . . 
Since this is the only parity violating term, we pick terms in the expansion of the log 
which have an odd power of Cq. Since we work to lowest order in Cq, this means we keep 
exactly one power, i.e. we simply renormalize the operator H^. Using Monge gauge we 
find that 

c * f r 

H * = ~2 / d2 4 W 12+^( W 22- W ll)+^ 7 ^l W l7 _W 2 W 2 7 )-I W 7 ( W l W 2 7 + W 2 W l 7 )- 2 ^ 2w 12j+--- 5 

(4) 

where u = d u = du/d^ 1 etc. and the ellipsis denotes terms with at least five fields or 
at least three 0's. While H eS is a complicated power series in u, 6, we can unambiguously 
determine c* ff by expanding i? eff to first order in u and zeroth order in 9, since the first 
term breaks parity and does not appear in any of the total derivative terms dropped in 
([^). While u 12 is a total derivative, so that this term seems to vanish, we can still compute 
it by giving the coupling Cq a small fictitious spatial dependence in intermediate stages of 
the calculation. 

We now quote the results of the calculation. Letting D = ^^log6 _1 , H^ eS looks 
like with Cq replaced by Cq ( 1 — ^ ) • Hence we find that the effect of fluctuations may 
be summarized by omitting them but replacing Cq by c* ff , where 

dct s (b-') _ - c ; s t 

dlogfr -1 ^^oTi 



(5) 



In general this H e s does not have the same functional form as H; in general we need to rescale 
the fields to recover the original relations between terms of different degrees of homogeneity (as 
for example in the O(n) nonlinear sigma model |ff?|| ). In our problem there is no need for an 
additional field renormalization of u. Intuitively, this is clear from the Monge gauge condition: 
a rescaling of a; 3 would by rotation invariance entail rescaling x 1 ' 2 which would spoil the gauge 
condition. Detailed calculation confirms that the nonlinear structure of (0)-(|3|) is retained if we 



do not rescale u [11|. 
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The constants k '7i were denned in (|l|). Strictly speaking they too should be allowed to 
run, but as mentioned they will arrive at the fixed line || («7i) e fj = 4K eff . Solving (|5|) we 
find 

<Ab- 1 ) = (b-T T/4 ™ e "c*o<c* , (6) 

the promised thermal softening. 

Eqn. (H) should not be construed as a temperature-dependence of c* ff , since in general 
the bare couplings will have some unknown T-dependence. However we can draw an 
interesting qualitative conclusion. The diameter R of helical ribbons (and possibly tubules 
as well) is controlled by a competition between chirality and stiffness. Oversimplifying 



somewhat by omitting the V x m boundary term (its effect is similar to that of ||12|j ), 
one finds in mean- field theory that |12| R oc Kq/cq. Since Cq can be much smaller than the 
cutoff A, R can be very large; indeed experimentally R can be ~ .5/xm ||. Hence thermal 
fluctuations can significantly modify the mean- field result; we can approximately account 
for their effects by writing R oc « eff /c* ff (Ail). Thus (|6|) says that varying Cq for fixed 
k , 7 1; T we have that the cylinder radius scales not as R oc (cq)" 1 but as 

Roc ( c *)-(i+r/4™ cff ) _ (7) 

The nice feature of (0) is that it may be possible experimentally to control Cg, without 
changing significantly the other bare parameters, simply by diluting the chiral amphiphiles 



with similar but achiral analogs [[19]]. Thus (0) is potentially a rather clean test of renor- 
malization effects in rigid chiral membranes. Whether the range of dilutions admitting 
helices or tubules will be great enough, and T/4%K eS can be made large enough, to test 
the scaling law ((7|) remains to be seen. Even if not, @ may still be applicable to rippled 



phases 12 



We have seen how the constituents of a membrane can express their chiral nature at 
long scales through the development of in-plane tilt order, and how thermal fluctuations 
can reduce the effective value of the bulk chiral coupling constant and in turn affect the 
shapes of self-assembled structures. Near the trivial, low temperature, fixed point the 
chiral term is relevant, but we have seen how thermal fluctuations reduce its effective 
dimension. This raises the possibility of a critical chiral membrane when this term becomes 
marginal. Unfortunately this will not happen at weak coupling, so we can say little about 
this intriguing possibility. 
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